home *** CD-ROM | disk | FTP | other *** search
/ C/C++ Users Group Library 1996 July / C-C++ Users Group Library July 1996.iso / vol_300 / 334_01 / controls.dem < prev    next >
Text File  |  1991-02-04  |  742b  |  30 lines
  1. #
  2. # warning:  this demo is SLOW on PCs without math coprocessors!
  3. #
  4. # From _Automatic_Control_Systems_, fourth ed., figure 6-14
  5. # transient response of a second-order system to a unit step input function
  6. #
  7. damp(t) = exp(-s*wn*t)/sqrt(1.0-s*s)
  8. per(t) = sin(wn*sqrt(1.0-s**2)*t - atan(-sqrt(1.0-s**2)/s))
  9. c(t) = 1-damp(t)*per(t)
  10. #
  11. #    wn is natural undamped frequency
  12. #    s is damping factor
  13. #
  14. wn = 1.0
  15. set xrange [0:13]
  16. set samples 50
  17. set dummy t
  18. #
  19. # plot c(t) for several different damping factors s
  20. #
  21. plot s=.1,c(t),s=.3,c(t),s=.5,c(t),s=.7,c(t),s=.9,c(t),s=1.0,c(t),s=1.5,c(t),s=2.0,c(t)
  22. pause -1 "Hit return to continue"
  23.  
  24. # undo what we have done
  25. set xrange [-10:10]
  26. set autoscale xy
  27. set samples 160
  28. set dummy x
  29.  
  30.